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Abstract 

A number  of  topics  in  finite  elasticity  theory  which  appear 
to  lend  themselves  to  further  development  are  briefly  discussed. 
These  include  (i)  the  effect  of  kinematic  constraints  which  are 
exactly,  or  approximately,  satisfied;  (iij  the  mechanics  of 
elastic  membranes;  (iii)  the  applicability  of  results  in  finite 
elasticity  theory  to  problems  involving  stress  relaxing  materials; 
(iv)  the  development  of  necessary  and  sufficient  conditions  for 
material  stability  of  isotropic  elastic  materials;  (v)  the  con- 
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1 . Introduct ion 

The  object  of  this  paper  is  to  present  some  research 
directions  in  finite  elasticity  theory,  which  have  not,  in 
the  opinion  of  the  author  been  fully  exploited.  The  choice 
of  topics  has  been  governed  by  discipline  oriented  rather 
than  problem  oriented  considerations.  That  is  to  say,  the 
choice  has  been  determined,  not  by  the  importance  of  the 
engineering  applications  which  may  stem  from  the  researches 
discussed,  but  rather  by  a desire  to  indicate  directions  in 
which  the  theory  of  finite  elasticity  can,  in  its  present 
state  of  development,  be  extended  in  a logical  fashion. 

With  these  objectives  in  mind,  the  basic  equations  of 
finite  elasticity  are  succinctly  presented  in  §2.  Then,  in 
§ § 3 and  5,  the  modifications  to  these  equations  which  result 
from  the  existence  of  kinematic  constraints  on  the  possible 
deformations  of  the  material  are  discussed.  These  constraints 
may  be  local  in  character  (§3),  or  global  (§5).  In  §6,  we 
discuss  the  manner  in  which  problems  could  be  solved  in  situa- 
tions in  which  the  kinematic  constraint  conditions  are  only 
approximately  satisfied. 

In  §4,  we  draw  attention  to  the  applicability,  with  minor 
modifications,  of  certain  results  of  finite  elasticity  theory 
to  problems  involving  viscoelastic  solids  with  fading  memory. 
In  §7  we  mention  briefly  some  recent  ideas  on  the  mechanics  of 
elastic  membranes,  due  mainly  to  Skalak  and  his  collaborators, 
which  would  appear  to  merit  further  exploitation. 
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Finally  in  § § 8 and  9,  we  describe  some  recent  results  of 
Sawyers  and  Rivlin  on  material  stability  of  an  incompressible 
isotropic  elastic  material  and  on  the  conditions  for  a bifurca- 
tion from  the  pure  homogeneous  deformation  of  a thick  plate  of 
such  a material  to  exist  when  it  is  subjected  to  a thrust. 
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2 . The  basic  equations  of  finite  elasticity  theory 

We  consider  the  isothermal  deformation  of  a body  of 
elastic  material.  Some  identifiable  configuration  of  the 
body,  generally,  but  not  necessarily,  that  in  which  the 
material  is  unstressed,  is  taken  as  the  reference  config- 
uration. We  identify  the  particles  of  the  body  by  their 
vector  positions  X , with  respect  to  a fixed  origin,  in 
the  reference  configuration.  As  the  body  is  deformed,  the 
vector  positions  of  its  particles  change  to  x , where  x 
depends  on  X and  the  time  t , thus: 

x = x (X , t)  . [2.1] 

% 

We  denote  the  components  of  X and  x in  a fixed  rectangular 
cartesian  coordinate  system  x by  XA  (A=l,2,3)  and  x.^  (i=l,2,3) 
respectively. 

The  strain-energy  per  unit  mass,  denoted  W , is  assumed 
to  depend  on  the  local  deformation  gradients  x.  . . For 
brevity,  we  will  use  the  notation 


g = 


giA 


= x 


i ,A 


[2.2] 


g is  called  the  deformation  gradient  tensor . 

Since  the  superposition  on  the  assumed  deformation  of  a 
rigid  rotation  leaves  W unchanged,  it  follows  from  rather 
simple  mathematical  considerations  that  W must  be  express- 
ible in  terms  of  the  deformation  gradients  through  the 
components  in  the  system  x of  the  symmetric  tensor  C , 

the  Cauchy  strain  tensor,  defined  by 
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~ = %+g~  * [2.3] 

where  the  dagger  denotes  the  transpose.  Thus, 

W = W(C)  . [2.4] 

Symmetry  of  the  material  in  its  reference  state  implies 
restrictions  on  the  form  of  W . Thus,  if  S denotes  a 
generic  transformation  of  the  3*3  matrix  representation  of 
the  symmetry  group  S of  the  material,  W must  satisfy  the 
relation 

W(SCS+)  = W(C)  . [2.5] 


The  group  S and  correspondingly  its  3x3  matrix  representation, 
is  either  the  full  orthogonal  group  or  a sub-group  of  it. 
Equation  [2.5]  must  be  satisfied  for  all  elements  of  the 
group  S . 


In  the  case  when  S is  either  the  full  or  proper  orthogonal 
group,  corresponding  respectively  to  an  isotropic  material  which 
is  or  is  not  centrosymmetric  in  the  reference  configuration,  it 
follows  that  W must  be  expressible  as  a function  of  I , I 
and  I ^ , where 


11  - tr  C , I2  = j[(tr  C)2  - tr  C2]  , 
1^  = det  C , 


[2.6] 


thus : 


W = WCI1,I2»I3)  • 


[2.7] 


For  deformations  which  are  possible  in  a real  material, 
det  g (which  physically  is  the  ratio  between  the  volumes  of 
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a material  element  at  time  t and  in  the  reference  configura- 
tion) must  be  positive  and  it  follows  that 

det  g - I*  . [2.8] 

Accordingly,  W is  expressible  as  a function  of  ^ , 1 , 
det  g , thus: 

W = W(I1,I2,det  g)  . [2.9] 


For  materials  which  are  transversely  isotropic,  the 
restrictions  implied  by  eq.  [2.5]  were  obtained  by  Ericksen 
and  Rivlin  (1).  Those  corresponding  to  the  cases  when  S 
is  one  or  other  of  the  crystallographic  point  groups  were 
determined  by  Smith  and  Rivlin  (2) . 


The  Cauchy  stress  tensor  a and  Piola-Kirchhoff  stress 
tensor  n are  given  in  terms  of  W by* 


~ = pI3^§f)  ’ 5 


f 3 W-j 

p(s?) 


t 


[2.10] 


where  p is  the  material  density  in  the  reference  configuration. 


a and  II  must,  of  course,  satisfy  appropriate  equations 
of  motion  which  are 


X3~x  * ~ + p$  = » Vn~  + p<£  = p*  » [2.11] 

where  denotes  the  spatial  gradient  taken  with  respect  to 

the  configuration  at  time  t and  V denotes  that  taken  with 

respect  to  the  reference  configuration. 

_ 

By  9W/3g  , we  mean  the  matrix  with  element  3W/3x. 
the  ith  row  and  Ath  column.  ’ 


in 
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The  stress  boundary  conditions  on  a and  II  are 
f ® on  , F = n^N  , 


[2.12] 


where  n and  N are  the  unit  normals  to  the  boundary  at 
time  t and  in  the  reference  configuration  respectively 
and  f and  F are  the  applied  forces,  per  unit  surface 
area  measured  at  time  t and  in  the  reference  configuration 
respectively. 

If  the  material  considered  is  isotropic,  so  that  W 
may  be  written  in  the  form  given  in  eq.  [2.7],  then  from 
eqs.  [2.3],  [2.6],  [2.7],  and  [2.10]  it  follows  that  a 
is  given  by 

a = 2Pl3Js[  (W1+I1W2)c  - W2c2  + ^36]  , [2.13] 


where  W.,  denotes  9W/9I  (a=l,2,3)  , 6 denotes  the  unit 

Ufc  u.  ~ 

tensor,  and  c is  the  Finger  strain  tensor  defined  by 


c 


[2.14] 


Apart  from  the  explicit  forms  for  W for  transversely 
isotropic  materials  and  for  materials  with  symmetries  corre- 
sponding to  one  or  other  of  the  crystallographic  groups,  the 
results  discussed  so  far  were  known  very  early  in  the  present 


century. 


If  0 (i=l,2,3)  are  the  principal  stresses  and  X.^ 

(i=l,2,3)  are  the  principal  extension  ratios,  then  I , Ig, 
I^  may  be  expressed  by 


X2  + X?  + X2 


= X 2X  2 + 
2 3 


X2X2 

A3A1 


X2X2 
1 2 


T ~ X2X2X2 
X3  ~ A1A2A3 


[2.15] 
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3,  Kinematic  constraints  and  controllable  deformations 

If  the  material  considered  is  incompressible,  then  the 
deformations  which  are  possible  in  it  are  subject  to  the 
constraint 

I3  = 1 [3.1] 

and  correspondingly  the  expression  for  a in  eq.  [2.10]  is 
replaced  by 

o = Pg(||)  " p6  , [3.2] 

where  p is  arbitrary  if  the  deformation  is  specified, 
reflecting  the  fact  that  the  superposition  of  an  arbitrary 
hydrostatic  pressure  on  the  force  system  applied  to  the  body 
leaves  the  deformation  unaltered.  In  the  case  when  the 
material  is  isotropic,  eq.  [2.8]  becomes 

W = W(I],I.o)  [3.3] 

and  from  eq.  [3.2]  we  find  that  eq.  [2.13]  is  replaced  by  (3) 

o = 2p[(W1+I1W2)c  - W2c2]  - p6  . [3.4] 

The  effect  of  constraints  on  the  expression  for  the  Cauchy 
stress  was  investigated  further  by  Ericksen  and  Rivlin  (1). 

They  considered  the  possibility  that  a number  of  constraints 
on  the  deformation  of  the  form 

fa(g)  = 0 (o-l,2,...)  [3.5] 

are  simultaneously  present  in  the  material  and  showed  that 
eq.  [2.10]  for  the  Cauchy  stress  must  then  be  replaced  by 
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where  the  p's  are  arbitrary  if  the  deformation  is  specified. 

An  example  of  such  constraints  is  provided  by  the  case 
when  there  are  one  or  more  directions  of  inextensibility  in 
the  undeformed  material.  If  elements  of  length  which  are 
initially  in  the  direction  of  the  x^  axis  are  inextensible , 
the  corresponding  constraint  condition  is 


The  case  when  there  are  two  directions  of  inextensibility  and 
in  addition  the  material  is  incompressible  was  first  studied  by 
Adkins  and  Rivlin  (4).  This  case  and  that  in  which  the  material 
is  reinforced  with  only  one  set  of  cords  was  further  studied  by 
them  and  by  a number  of  other  workers,  notably  by  Spencer,  Pipkin 
and  Rogers.  Excellent  accounts  of  much  of  this  work  is  given  in 
the  monograph  by  Spencer  (5)  and  the  review  article  by  Pipkin  (6). 

The  constraint  imposed  by  incompressibility  is,  as  is  now 
well-known,  much  more  significant  in  its  implications  than 
might  at  first  sight  appear.  For,  it  enables  us  to  solve 
certain  problems  without  making  any  specific  assumptions 
regarding  the  strain- energy  function  than  those  implied  by 
isotropy  and  incompressibility  and  expressed  by  eqs.  [3.3] 
and  [3.1].  Although  these  problems  are  admittedly  simple  ones 
they  lend  themselves  to  experimental  realization  and  enable 
one  to  take  the  critical  step  of  actually  determining  for  a 
particular  material  - vulcanized  rubber  - the  manner  in  which 
the  strain-energy  function  depends  on  1 ^ and  1^  and  to  verify 


the  predictive  value  of  the  strain- energy  function  so  determined. 
For  example,  if  we  consider  a thin  square  sheet  with 


thickness  h and  unit  edge,  it  is  fairly  easy  to  show  that 


the  forces  T 

and  T2  , 

which  must  be  applied  in 

order  to 

deform  it  into 

a rectangle  of  dimensions  A and 

A2  , are 

given  by  (7) 

Tx  = 2ph| 

1 2 

[3.8] 

T2  = 2ph| 

(X2  - a31x2)(W1*X1W2)  ' 

where 


- a2  ,2 

A1  A2 


X2X2 

A1A2 


i = — + — + a2a2 

2 ,2  ,2  12 

A1  A2 


[3.9] 


Experimentally,  we  can  measure  T and  T0  for  specified 

values  of  A and  A and  then  use  eqs.  [3.8]  to  determine 
12 

W1  and  W2  for  the  corresponding  values  of  I1  and  I2 
obtained  from  eqs.  [3.9]. 


Again,  we 
of  radius  a , 
twist  per  unit 
the  torsion  we 
forces  0 per 
given  by 


may  consider  the  torsion  of  a circular  cylinder, 
by  forces  applied  to  its  plane  ends.  If  the 
length  is  ip  , we  find  that  in  order  to  maintain 
have  to  apply  to  the  plane  ends  circumferential 
unit  area  and  normal  forces  Z per  unit  area, 


Z 


0 = 2ipr p (W1  + W2)  , 

-2ip2pfr2W2  + f rW^dr 

' ' a 


[3.10] 


where,  in  this  case, 


l] 


I 


l1  = I2  = 3 + 4>2r2  . [3.11] 

Deformations  which  can  be  supported  without  application 
of  body  forces  and  for  which  the  surface  forces  can  be  calcu- 
lated without  making  specific  assumptions  regarding  the  form 
of  the  strain-energy  function,  beyond  those  implied  by  such 
material  symmetry  and  kinematic  constraints  as  may  be  relevant 
to  the  class  of  materials  under  consideration,  are  called 
controllable  deformations . The  attempt  to  find  all  possible 
controllable  deformations  in  an  isotropic  material  was  initiated 
by  Ericksen  (8)  and  considerably  extended  by  Singh  and  Pipkin  (9) 
and  by  Klingbeil  and  Shield  (10). 

In  the  case  of  compressible  isotropic  elastic  materials, 
it  was  shown  by  Ericksen  (11)  that  the  class  of  controllable 
deformations  consists  only  of  homogeneous  deformations.  In 
the  case  of  incompressible  materials  it  is  substantially  wider. 

We  may  also  identify  a class  of  deformations  which  may  be 
called  quasi -controllable  deformat  ions . As  an  example  of  these, 
we  may  consider  (12)  an  incompressible  isotropic  elastic  material 
to  be  contained  in  the  annular  region  between  two  rigid  coaxial 
circular  cylinders,  one  of  which  is  moved  parallel  to  the  other 
by  a fixed  amount,  or  rotated  relative  to  the  other,  about  the 
common  axis,  through  a specified  angle.  In  either  case  the 
displacements  of  the  particles  of  the  material  are  determined 
as  a function  of  radial  position  r by  a single  ordinary 
differential  equation  which  depends  on  the  form  of  W . The 
surface  forces  which  must  be  applied  in  order  to  maintain  the 
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deformation  can  then  be  calculated  leaving  open  the  precise 
form  of  the  strain-energy  function  and  the  consequent  radial 
distribution  of  the  displacement. 

Another  type  of  system  which  leads  to  tractable  problems 
is  the  elastic  membrane  - typically  a vulcanized  rubber  membrane  - 
which  can  be  subjected  to  large  deformations.  Work  on  such 
systems  was  first  carried  out  by  Adkins  and  Rivlin  (13).  If 
the  membrane  and  the  deforming  forces  have  rotational  symmetry 
about  an  axis,  the  local  deformation  is  a pure  homogeneous 
deformation  with  principal  directions  along  the  lines  of  longi- 


tude and  latitude  in  the  deformed  membrane  and  principal 


extension  ratios  in  the  surface  of  the  membrane,  X , X^  say. 
The  corresponding  two-dimensional  Piola-Kirchhoff  stress  is 
then  given  by  eqs.  [3.8],  If  W is  any  known  function  of  1^ 
and  I2  , the  point  equations  of  equilibrium  can  be  integrated 
numerically. 
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4 . An  application  of  finite  elasticity  theory  to  inelastic 
materials 

It  was  noted  by  Rivlin  (14)  that  the  expressions  obtained 
for  the  forces  which  must  be  applied  to  maintain  controllable 
deformations  in  isotropic*  elastic  materials  may  be  applied, 
with  changes  which  are,  from  a mathematical  standpoint,  trivial 
in  nature,  to  obtain  expressions  for  the  forces  which  must  be 
applied  in  order  to  maintain  identical  deformations  in  stress- 
relaxing  materials. 

We  consider  that  a body  of  stress -relaxing  material  is 
deformed  during  a small  time  interval  AtQ  , say,  at  the 

instant  of  time  t , say,  and  that  the  deformation  is  then 

o 

held  constant.  We  assume  that  after  the  lapse  of  a sufficient 
time,  the  Cauchy  stress  a depends  on  the  deformation  gradient 
g and  the  time  that  has  elapsed**,  t say,  after  the  deforma- 
tion is  produced.  Thus, 

o = o(g,t)  . [4.1] 

It  can  easily  be  shown,  from  the  consideration  that  if 
an  arbitrary  rigid  rotation  is  superposed  on  the  assumed 
deformation,  the  Cauchy  stress  tensor  is  rotated  by  an  equal 
amount,  that  a must  be  expressible  in  the  form 

a = gf  (C,t)g+  , [4.2] 

Similar  considerations  apply,  in  principle,  to  stress-relaxing 
anisotropic  materials. 

** 

The  extent  to  which  this  assumption  is  valid  depends,  of  course, 
on  the  nature  of  the  st re s s -relax i ng  material  and  on  the  class 
of  allowable  deformation  histories  by  which  the  final  state 
of  static  deformation  is  attained. 
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where  C is  the  Cauchy  strain  defined  by  eq.  [2.3].  This 
step  in  our  argument  parallels  that  by  which  the  passage  was 
achieved  from  the  assumption  for  an  elastic  material  that  W 
is  a function  of  g to  the  result  expressed  by  eq.  [2.4]. 
Again,  paralleling  the  passage  from  eq.  [2.4]  to  eq.  [2.5], 
we  find  that  if  the  material  has  symmetry  in  its  reference 
state,  then  with  the  notation  of  §2 

f(SCS+)  = Sf(C)S+  . [4.3] 


For  an  isotropic  material  this  leads  to  the  conclusion  that 
o must  be  expressible  in  the  form 

p 

a = a^c  + a2c  + aQ6  , [4.4] 

where  the  a's  are  functions  of  1^,  I2  and  1^  and  of  t and 
c is  the  Finger  strain  defined  by  eq.  [2.14].  For  an  incom- 
pressible material,  this  relation  becomes 


2 

a = a c + a2c  - p5  , [4.5] 

where  the  a's  are  now  functions  of  I , I0  and  t and  p 
is  an  arbitrary  hydrostatic  pressure.  Comparing  this  equation 
with  eq.  [3.4],  we  see  that  the  latter  can  be  obtained  from  eq. 
[4.5]  by  replacing  a^  and  a ^ by  2p(Wi+I^W2)  and  -2pW2 
respectively. 

The  forces  necessary  to  maintain  a controllable  deformation 
is  a stress-relaxing  material  can  be  obtained  from  those  for  an 
elastic  material,  by  replacing  W and  W0  by  (a1+I1a2)/2p 
and  -a2/2p  respectively.  For  example,  the  formulae  [3.10]  for 
the  azimuthal  and  normal  forces,  0 and  Z per  unit  area,  which 
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must  be  applied  to  the  plane  ends  of  a right -circular  cylinder 
in  order  to  maintain  in  it  a simple  torsion  of  amount  ip  become 


0 = i^r  [ot1  + [2  + iJ;2r2)a2] 


-a, 

C 


r2  + 


r{a1  + (3  + i^r"  )a?  }dr 


[4.6] 


2 2 


where  the  a's  are  functions  of  ip  r and  of  elapsed  time 


after  the  deformation  is  produced. 
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5 . Global  constraints 

A somewhat  different  type  of  constraint  from  those  discussed 
in  §3  may  be  appropriate  to  certain  systems.  For  example,  if  a 
body  of  vulcanized  rubber  is  swollen  homogeneously  in  some  solvent 
and  is  then  inhomogeneously  deformed,  the  solvent  may  migrate 
from  one  part  of  the  body  to  another.  Since  both  the  vulcanized 
rubber  and  the  solvent  separately  may,  with  good  approximation, 
be  treated  as  incompressible,  the  total  volume  of  the  body  will 
be  substantially  unaltered  by  the  deformation,  but  each  element 
of  the  vulcanized  rubber  may  change  the  volume  it  occupies  due 
to  the  absorption  or  loss  of  solvent.  Accordingly,  1^  , defined 
in  eqs.  [2.6],  may  be  regarded  as  a measure  of  the  volume  frac- 
tion of  solvent  at  any  point  and  the  strain-energy,  at  any  rate 
for  quasi-static  deformations  of  the  swollen  rubber  vulcanizate 
must  be  expressible  in  the  form  [2.7].  The  constraint  implied 
by  the  constancy  of  volume  of  the  body  may,  of  course,  be  written 


det  g • dV  = V , 


[5.1] 


where  V is  the  domain  occupied  by  the  body  in  the  homogeneous 
reference  configuration.  The  point  equations  which  can  be 
derived  for  the  type  of  system  considered  are  similar  to  those 
for  a compressible  material,  except  that  the  force  boundary 
condition  is  now  given  by 


[5.2] 


and  the  constraint  condition  by  eq . [5.1]. 
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6.  Perturbations  of  the  constraint  condition 


We  note  that  if  we  superpose  on  the  deformation  of  a body 
described  by  eq.  [2.1],  a uniform  dilatation,  so  that  x is 
replaced  by  Xx  , where  X is  independent  of  X , then  the 
Finger  strain  c defined  by  eq.  [2.14]  is  replaced  by  X c . 


Similarly,  the  strain  invariants  1^  and  , defined  by 

P h f) 

eqs.  [2.6],  are  replaced  by  X 1^  , X I0  and  XI  respectively, 
Thus,  c,  and  J2  defined  by 


c = i * 1/3  c J = i i ' 1//3 
2 3 ~ ’ J1  il13 


J.  = I„I. 


[6.1] 


r 3 ~ * 1 13  ’ 2 2 3 1 J 

are  unaltered  by  the  superposed  dilatation.  This  suggests 
(Flory  (15),  Penn  (16))  that,  for  an  isotropic  material,  we 
express  W as  a function  of  and  , thus 

W = W(J1,J2,I3)  , [6.2] 

rather  than  as  a function  of  1.^,  I2,  as  indicated  in  eq. 

[2.7].  More  conveniently  for  our  purposes,  we  can  express  W 
as  a function  of  J , J2  and  x (=1^)  , thus: 


W = W(J1,J2,x)  . 

Then,  the  expression  [2.13]  for  a can  be  rewritten  as 

2p  17  9W  , T 3W  \-  9 W ~2 

2 = + J1  372Js  ' 9T2  2 

■ 1 T 3W  2 T aw  . 1 T 3W  1 A 

['7  jl  ' I 3 2 3J2  1 T 97  J • 


[6.3] 


[6.4] 
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W = wQ(JlfJ2)  + Cx-l)w1(J1,J2)  + \ (t-1)2w2(J1,J2)  , [6.5] 


where 


WU  = 1 * w1(J1>J2) 


aw 

3r 


T = 1 


f t T , 32W 

W2(J1’J2>  = 77 


[6.6] 


T=1 


Then,  if  we  introduce  eo.  [6.5]  into  eq.  [6.4]  and  neglect 
terms  of  higher  degree  than  the  first  in  t-1  , we  obtain 
3w„  3w 3w, 


= - *R  * ^ - 


3 J 


0 ^2 
c 


1 , 3w0  + 2 . 9"o  1 .. 

i Ji  7 I J2  w;  ' i wi|5 


3 w. 


3 w 


. 3w  3 w 

+ 2p(T'1)  {{  37  + J1  3R  ' U7  ' J1  3J2 


°'a 


/Sw  0 


( 1 0V2  [1  . 

' 1^7  ' 57; j ; [j  - 5J1 


2 T 

fdwo 

3w  s 

1 

+ 3 J2I 

k ' 

3Rj 

+ 2 W2 

7 

[6.7] 

We  consider  a sequence  of  constitutive  equations  for  nearly 
compressible  materials,  with  specified  density  p in  the  undeformed 
state,  which  has  as  its  limit  the  constitutive  equation  for  an 
incompressible  material  with  strain-energy  function  given  by 

W = w0CJlfJ2)  . [6.8] 

The  corresponding  expression  for  the  stress  is 

3w, 


; • S)- 


0 -21 
377  ? 1 


p6  , 


[6.9] 
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where  p is  an  arbitrary  hydrostatic  pressure.  We  note  that 
since  all  deformations  in  an  incompressible  material  are 
necessarily  isochoric,  we  have  c=c  in  this  case. 

The  fact  that  for  the  incompressible  material  the  con- 
stitutive equation  for  o must  take  the  form  [6.9]  leads  to 
the  restriction  on  the  coefficient  of  2p(x-l)  in  eq.  [6.7] 
that  it  tends  to  infinity,  in  such  a way  that  its  product  with 
(t-1)  remains  finite,  as  the  material  tends  to  incompressibility. 

We  can  solve  problems  for  slightly  compressible  materials, 
for  which  eqs.  [6.5]  and  [6.7]  are  valid,  by  first  solving  the 
problem  for  an  incompressible  material  for  which  the  eqs.  [6.8] 
and  [6.9]  are  valid.  Then,  using  a perturbation  procedure,  we 
can  find  the  correction  to  the  displacement  field  appropriate 
to  the  slightly  compressible  material. 

Analogous  procedures  should  be  possible  in  situations  where 
other  constraint  conditions  than  incompressibility  are  approxi- 
mately valid.  However,  it  must  be  borne  in  mind  that  the 
perturbations  involved  need  not  necessarily  be  regular  (Rogers 
and  Pipkin  (17),  Everstine  and  Pipkin  (18)). 


; 
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7.  Two-dimensional  materials 


One  can  construct  an  interesting  two-dimensional  analog 
to  finite  elasticity  theory,  which  is  apparently  applicable 
to  various  living  tissues,  e.g.  the  membrane  of  the  human  red 
blood  cell  (Hochmuth  and  Mohandas  (19),  Skalak,  Tozeren,  Zarda 
and  Chien  (20) ) . 

We  consider  a membrane  of  transversely  isotropic  elastic 
material,  the  axis  of  rotational  symmetry  being  normal  to  the 
membrane.  With  the  usual  assumptions  of  membrane  theory,  we 
can  write  the  strain-energy  per  unit  area  (measured  in  the 
reference  state)  as  a function  of  I , 10  defined  in  terms  of 
the  principal  extension  ratios  A^,  Ao  by 


? ? ? o 

h ’ A'l  * ‘s  ' 's  ' Vi  ' 


(7.1] 


thus : 


w = w(i1 ,i2; 


[7.2] 


The  principal  line  stresses  a (a=l,2)  arc  then  given  by 


a =21 

a 


.2  5W  . aw 

a 51  l2  rro 


[7.5] 


It  appears  that  for  certain  tissues,  the  area  of  an  element 
of  the  membrane  reamins  approximately  constant  as  the  membrane 
is  deformed.  We  are  therefore  led  to  construct  a constitutive 
equation  which  idealizes  this  observation  by  imposing  the 
constraint  condition* 


I2  = 1 [7.41 

If  the  material  of  the  membrane  If  also  incompressible,  this 
implies,  of  course,  that  it  is  inextensible  in  a direction 
normal  to  it . 
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on 


the  deformation.  Then,  W is  a function  of  1^  only,  thus: 

W = W(I  ) [7.5] 

and  the  constitutive  equation  [7.3]  is  replaced  by 


3 W 

0 —2  a.  *ry  “ P 

a a 9L 


[7.6] 


where  p is  arbitrary  if  the  deformation  is  specified. 

If  the  constraint  [7.4]  is  only  approximate,  the  resulting 
change  in  the  deformation  could  be  calculated  in  a manner  analogous 
to  that  introduced  in  the  case  of  slightly  compressible  materials. 


; 1 
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8 . Initial  stress  problems 

We  have  seen  that  for  incompressible  isotropic  elastic 
materials,  certain  simple  problems  can  be  solved  without 
particularizing  the  form  of  the  strain-energy  function  beyond 
that  given  in  eq.  [3.3].  The  recognition  of  this  fact  led  to 
the  formulation  of  a rational  theory  for  the  solution  of  initial 
stress  problems  in  which  the  strain-energy  function  is  left 
in  the  form  [3.3]  and  the  initial  stress  is  that  corresponding 
to  one  or  other  of  the  so-called  controllable  solutions. 

Insofar  as  controllable  solutions  exist  for  compressible 
isotropic  materials,  analogous  considerations  apply. 


The  Piola-Kirchhoff  stress  IT  for  a compressible  material 
corresponding  to  a deformation  X+x  is  given  by  eq . [2.10], 

viz: 


n 


Ai 


3W 

9Xi,A 


[8.1] 


We  now  consider  an  infinitesimal  superposed  deformation 
x -*■  x+eu  , where  e is  sufficiently  small  so  that  we  can 
linearize  in  it.  This  superposed  deformation  may  be  time- 
independent  or  time -dependent . The  incremental  stress  eiTAi 
is  then  given  by 

*2„ 

u.  „ . [8.2] 


a2w 

^Ai  P dx.  ,3x. 


J.B 


i > A j ,B 

If  the  material  is  isotropic,  we  can  introduce  the  expression 
[2.7]  for  W and  carry  out  the  indicated  differentiations  to 
obtain  the  appropriate  expression  for  . In  this  expression 

we  substitute  for  the  coefficients  of  u , the  expressions 

J 


J 


! 


I 


'■! 

( 

i i 


u 


'■f 


H 

* 
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appropriate  to  the  known  underlying  static  deformation.  The 
appropriate  equation  of  motion  for  the  superposed  deformation 
is 


%i,A  = pUi 


[8.3] 


and  the  force  boundary  condition  is 


TT  . . N 
Ai  A 


K > 


[8.4] 


where  ef.  is  the  additional  force  associated  with  the  super- 
posed deformation. 


In  the  case  when  the  material  is  incompressible,  the 
expression  [8.2]  for  the  incremental  stress  must  be  modified 
to  reflect  the  fact  that  the  superposition  of  a hydrostatic 
pressure  does  not  change  the  deformation.  In  addition,  we 
have  to  add  an  equation  which  expresses  the  fact  that  no  volume 
changes  can  take  place. 

On  the  basis  of  the  considerations  outlined  above,  a 
number  of  initial  stress  problems  have  been  solved.  The 
earliest  of  these  was  the  calculation  by  Green  and  Shield  (21) 
of  the  effect  of  an  initial  finite  extension  on  the  torsional 
modulus,  for  infinitesimal  torsion,  of  a rod  of  arbitrary  but 
uniform  cross-section.  This  generalized  an  earlier  result  of 
Rivlin  (12)  for  a rod  of  circular  cross-section.  The  result 
of  Green  and  Shield  (21)  was  corrected  and  verified  experi- 
mentally by  Gent  and  Rivlin  (22).  Shortly  afterwards,  Green, 
Rivlin  and  Shield  (23)  discussed  the  effect  of  initial  stress 
on  the  indentation  by  a punch  of  an  elastic  half-space  subjected 
to  an  initial  biaxial  pure  homogeneous  deformation  with  two 
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equal  principal  extension  ratios  which  define  a principal  plane 
parallel  to  the  bounding  plane  of  the  half-space.  These  results 
and  many  others  of  a similar  character  have  been  discussed 
extensively  in  the  subsequent  literature.  Here  we  shall  con- 
sider only  the  initial  stress  problem  for  plane  waves  of 
infinitesimal  amplitude  propagating  in  a deformed  elastic 
materi al . 

First,  let  us  suppose  that  the  medium  in  which  the  plane 
wave  propagates  is  a space  of  isotropic  elastic  material  sub- 
jected to  a pure  homogeneous  deformation  with  principal 
extension  ratios  A , X^,  and  principal  directions  parallel 

to  the  axes  of  a rectangular  cartesian  coordinate  system  x . 

Let  us  suppose  also  that  the  wave  is  sinusoidal,  has  angular 
frequency  w , wave  number  k , slowness  S = k/w,  and  prop- 
agates in  the  direction  of  the  unit  vector  n . Then,  with 
the  usual  complex  notation,  we  can  express  the  displacement 
vector  eu  associated  with  this  wave  in  the  form 

u = A exp  t(kn-x  - ut)  . [8.5] 

The  secular  equation  determining  the  complex  slowness  S 

was  obtained  by  Hayes  and  Rivlin  (24)  in  the  case  when  the 

material  is  compressible.  Here,  however,  we  shall  give  the 

, somewhat  simpler  results,  obtained  by  Sawyers  and  Rivlin  (25), 

, 

for  the  case  when  the  material  is  incompressible.  In  this  case 


the  secular  equation  is 
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a=  4(X?n^+  — ) { (W  + A*W  _)  (W  + X*W  _)  + ...] 


, 2_ 2 


11 


+ Wi[n2n3^A2‘A3)2(Wl+AlW2)Al  + •••!) 

+ 4W2[X^X^n2n2(X2-X3)2(Wl+x2w2)Al  + ...] 

- 16n2n2n2 (X2  - X2 ) 2 (X2  - X2) 2 (X2  - X2)  2 (W22  "W^W^g) 


B = 2(W1+x2w2)  [X2n2  + X2n2  * (X2 -X^n^n^ ; 


[8.7] 


where  the  dots  represent  terms  obtained  from  those  given  by 
cyclic  permutation  of  the  subscripts  on  the  X's,  n's  and 
A' s and  where  the  A's  are  defined  by 


A,  = 2(X, 


V 


wx+x2w2 


[8.8] 


with  analogous  expressions  for  A0  and  AQ  . WX1  > w12  and 
W22  denote  92W/9I2,  92W/3I  31  and  32W/9I2  respectively. 

2 

We  may  regard  eq.  [8.6]  as  a quadratic  equation  in  S 

The  necessary  and  sufficient  conditions  for  it  to  yield  two 

. 2 
positive  roots  for  S are 


a>0  , B >0  , B2  > 4a  . 


[8.9] 


Following  Hadamard, we  adopt  these  conditions  as  the  conditions 
for  material  stability. 


It  would  be  of  interest  to  obtain  conditions  on  W such 
that  the  relations  [8.9]  are  satisfied  for  all  n , i.e.  for 
all  directions  of  the  normal  to  the  wave-front.  So  far  this 
has  not  been  achieved.  However,  Sawyers  and  Rivlin  (25)  have 
obtained  necessary  and  sufficient  conditions  on  W , so  that 
the  relations  [8.9]  are  valid  for  all  waves  for  which  n lies 
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in  a principal  plane  of  the  pure  homogeneous  deformation. 
In  this  case  necessary  and  sufficient  conditions  that  the 
relations  [ 8 . 9 j be  satisfied  are 


W 


11 


2XiWl? 


W1  + XiW2  > 


x>22 


t(!i  -4  -l) 


[8. io ; 


(i=l,2,3)  . 


The  conditions  [8.10]  were  previously  obtained  by  Baker  and 
Ericksen  (26)  and  are  usually  referred  to  in  the  literature  as 
the  Baker-Ericksen  conditions.  It  has  also  been  shown  by  Sawyers 
and  Rivlin,  in  unpublished  work,  that  the  condition  B > 0 
follows  from  the  Baker-Ericksen  conditions. 


The  conditions  [8.10]  can  also  be  obtained  (27)  from  the 
following  considerations.  We  consider  the  material  to  be  sub- 
jected to  a pure  homogeneous  deformation  as  before.  We  now 
superpose  an  infinitesimal  simple  shear  for  which  the  plane  of 
shear  is  a principal  plane  of  the  pure  homogeneous  deformation 
and  the  direction  of  shear  is  that  of  a unit  vector,  K say, 
lying  in  this  plane.  We  consider  the  change,  as  a result  of 
this  superposed  simple  shear,  in  the  traction  in  the  direction 
of  K acting  on  unit  area  of  a plane  parallel  to  K and  to 
the  normal  to  the  principal  plane  considered.  The  ratio  between 
this  traction  and  the  amount  of  shear  is  the  incremental  shear 
modulus.  The  relations  [8.10]  are  necessary  and  sufficient 
conditions  for  this  to  be  positive  for  all  choices  of  K in 
a principal  plane. 
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9 . Bifurcation  solutions  for  a thick  plate  under  thrust 

An  example  of  another  type  of  initial  stress  problem  which 
can  be  solved,  without  imposing  restrictions  on  the  form  of  the 
strain-energy  function  beyond  those  implied  by  isotropy  and 
incompressibility,  is  provided  by  the  calculation  of  the  bifur- 
cation condition  for  a thick  rectangular  plate  of  isotropic 
incompressible  elastic  material,  subjected  to  biaxial  loading. 

Suppose  the  edges  of  the  plate  are  parallel  to  the  axes 
of  a rectangular  cartesian  coordinate  system  x and  that  the 
lengths  of  these  edges  are  2 , 2 21  ^ (28^  >>  2 8^,  . 

The  plate  is  subjected  to  a pure  homogeneous  deformation,  with 
extension  ratios  A 1 , A^  and  principal  directions  parallel 
to  the  axes  of  the  system  x , and  is  held  so  that  A^  remains 
fixed,  the  faces  perpendicular  to  the  2-direction  being  force- 
free.  The  values  of  A , A 2 are  then  varied  (in  accordance, 
of  course,  with  the  incompressibility  condition  A^2  = ^3^ 
and  the  critical  value  of  A (=  A^/A  ) is  determined  for  which 
a bifurcation  of  the  flexural  or  barreling  type  in  the  12-plane 
can  exist.  The  faces  initially  perpendicular  to  the  1-direction 
are  constrained,  so  that  they  remain  perpendicular  to  the  1 -direction 
and  the  tangential  surface  traction  is  zero  over  them.  While 
these  constraints  are  admittedly  somewhat  artificial,  they  have 
the  merit  that,  with  them,  the  problem  can  be  solved  exactly. 

We  denote  the  wave-length  of  the  superposed  flexural  or 
barreling  deformation  by  2ir/ft  , so  that 

ft  = mr/28..^  (n=l,2,...)  . 


pw  >i» 


[9.1] 
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Then,  the  critical  value  of  X , at  which  a bifurcation  from 

the  state  of  pure  homogeneous  deformation  becomes  possible,  is 

given,  for  specified  values  of  n (i.e.  of  the  mode  of  the 

superposed  flexure  or  barreling)  and  of  ’ by 

tanh  /ft2+X2ft2  \ fi 

tanh  ft (fi^  + x2n2  J ^ 

or  [9.2] 

tanh  nia2  ftt22  + X2Sl2\  n 

tanh  V2  ' ^ ' 

accordingly  as  the  superposed  deformation  is  of  the  flexural 
or  barreling  type.  fi  and  ^2  are  defined  by 

a2,n2  = fi2{  [1+ A2  + AC1- X)  2 ] ± [{1+X2  + A(1-X)2}2-4X2]J5}  , [9.3] 


where  A is  defined  by 


2(X1+X2): 


VX3W2 


'Wll*2x23Wl2tX3W22) 


[9.4] 


Equations  [9.2]  were  apparently  first  derived  by  Wesolowski 
(28). 

It  follows  from  eq.  [9.3]  that  formally,  at  any  rate,  we 
may  discuss  the  implications  of  eqs.  [9.2]  in  three  distinct 
cases : 


(i) 

and 

are 

both  positive. 

(ii) 

«2 

and 

are 

complex  conjugates, 

(iii) 

and 

fl2 

are 

both  negative. 

It  is  apparent  from  eq . [9.3]  that  case  (iii)  arises  if  and 
only  if 


r 
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* < -(^r)2  • 


case  (ii)  arises  if  and  only  if 


r A + 1 2 . . . , 

■(rnv  i a < -i  , 


[9.5] 


[9.6] 


and  case  (i)  arises  if  and  only  if 


A > -1  . 


[9.7] 


The  necessary  conditions  [8.10]  for  material  stability 
derived  earlier  imply  that  we  need  consider  only  cases  (i) 
and  (ii).  We  shall,  in  fact,  consider  only  case  (i). 

In  order  to  discuss  the  qualitative  implications  of 
eqs.  [9.2],  Sawyers  and  Rivlin  (29)  found  it  convenient  to 
rewrite  them  in  the  form 


sinh[  (2\^cosh  6)n]  _ + [ (1  + A c ) 2 + 4 sinh2tS]cosh  6 


sinh[ (2A^sinh  6)n] 


[-(1-A)  +4  cosh'<5]sinh  <5 


[9.8] 


where 


n = 


nvl, 

V = nr 


and  sinh26  = (1+A)(A-1)2  [9.9] 


Ii 


V 

k 


and  the  plus  and  minus  signs  in  eqs.  [9.8],  apply  to  super- 
posed deformations  of  the  flexural  and  barreling  types  respectively. 

Equations  [9.8]  may  be  regarded  as  equations  for  the 
determination  of  > and  hence  of  n^g/^x  > corresponding 

to  a specified  critical  value  of  A . Then,  provided  A > "1  » 


car, 
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i.e.,  6 is  real,  it  is  fairly  easy  to  prove  the  following 

theorems : 

(a)  If  A<U  , corresponding  to  the  situation  when  the  applied 
force  in  the  1-direction  is  tensile  or  zero,  neither  of 
the  eqs.  [9.8]  has  a real  solution  for  nl^/Z^  . Accord- 
ingly in  this  case  no  bifurcation  of  either  the  flexural 
or  barreling  type  is  possible. 

(b)  If  A>1  , corresponding  to  the  situation  when  the  applied 

force  in  the  1-direction  is  a thrust,  and  the  upper  sign 
in  eqs.  [9.8]  is  taken,  so  that  we  are  considering  only 
superposed  deformations  of  the  flexural  type,  then  one  and 
only  one  solution  for  exists  provided  that 

4 cosh2<$  - (A - 1 ) 2 > 0 [9.10] 

and  no  solution  exists  if 

4 cosh2 6 - (A  - 1 ) 2 < 0 . [9.11] 

(c)  If  A>1  and  the  lower  sign  is  taken  in  eqs.  [9.8],  so  that 
we  are  considering  superposed  deformations  of  the  barreling 
type,  then  one  and  only  one  solution  for  nZ^/Z^  exists 
provided  that 

4 cosh2  6 - (A  - 1 ) 2 < 0 [9.12] 


and  no  solution  exists  if 


4 cosh26  - (A - 1 ) 2 > 0 . [9.13] 

It  follows  that  whatever  may  be  the  strain-energy  function, 
subject  to  the  restriction  A>-1  , the  ranges  of  critical  values 
of  A at  which  bifurcations  of  the  flexural  and  barreling  types 
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can  occur  are  discrete.  The  value  of  A separating  these 
regions  is,  of  course,  given  by 

4 cosh26  - (A-l)2  = 0 , [9.14] 

and,  with  eq.  [ 9 . 9 ] 2 , by 

A3  - (3+A)A2  + (2A-1)A  - (1+A)  = 0 . [9.15] 

We  bear  in  mind  here  that  if  the  expression  for  W as  a function 
of  I1  and  I2  is  known  then  the  dependence  of  A on  A is 
also  known. 

For  any  constant  value  of  A we  can  plot  n against  the 
critical  value  of  A for  a bifurcation  to  occur.  We  obtain 
curves  of  the  form  shown  in  Fig.  1.  Generally,  the  situation 
A = constant  is  somewhat  artifical*.  However,  curves  such 
as  those  in  Fig.  1 can  be  used  to  determine  the  critical  values 
for  bifurcations  of  the  flexural  and  barreling  types,  corre- 
sponding to  any  specified  strain-energy  function  and  a specified 
value  of  n , and  hence  of  nl  /SL  . From  Fig.  1 we  read  off, 
for  the  specified  value  of  n » corresponding  values  of  A and 
A , for  flexure,  say.  We  plot  these  on  a graph  and  on  the  same 
graph  we  plot  corresponding  values  of  A and  A obtained  from 
the  expression  [9.4]  for  A with  the  specified  strain-energy 
function.  The  intersection  of  these  two  curves  gives  us  the 
critical  value  of  A for  a bifurcation  of  the  flexural  type 
to  exist.  An  analogous  procedure  using  A vs.  A curves  for 
barreling,  obtained  from  Fig.  1,  enables  us  to  determine  the 
_ 

An  exception  arises  in  the  case  of  the  Mooney-R ivl in  strain- 
energy  function  for  which  A=0  . 
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critical  value  of  X for  a bifurcation  of  the  barreling  type 
to  occur. 

A particularly  interesting  conclusion  results  in  the  case 
when  n , and  hence  n£2/^1  , is  small  (30).  Using  eq.  [9.8] 
with  the  plus  sign,  we  can  obtain  the  following  formula  for  the 
critical  value  of  n , corresponding  to  bifurcation  of  the 
flexural  type,  as  a function  of  n , in  the  form 

X-l  = a1n2  + a2ni*  + a3n8  + a^n8  +0(n10)  , [9.16] 

where 

ai  = I ’ a2  = re  ’ a3  = Tf  (A0  + ’ 

aU  “ 77  (A0  + TT?)  * 1 9 . 1 7 ] 

and 

A = Lt  A . [9.18] 

0 X-l 

It  ea  y t>  show  that  whatever  may  be  the  expression 
for  W , 

♦ M - 1 ) ‘ = A . ♦ 0 ( r/  ) . [9.19] 

We  see  from  I .it  that  the  critical  value  of  X-l  is 

I4  . 

independent  oi  \ up  to  terms  of  order  n . Also,  it  depends 
on  W only  through  it-  second  derivative,  evaluated  in  the 
state  for  which  the  thrust  i>  zero,  up  to  terms  of  order  n » 
so  that  for  all  materials  having  a strain-energy  of  the  Mooney- 
Rivlin  form,  the  critical  value  of  X-l  is  the  same  up  to  terms 
of  order  n8  . 

To  give  some  quantitative  feel  for  the  insentivity  of  the 
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critical  value  of  A to  the  value  of  n and  hence  of  * 

we  mention  that  if  AQ  = 0 and  n = .793  , which  for  n=l  corre- 
sponds to  = 0.505  , the  value  of  A obtained  by  taking 

the  first  three  terms  in  eq.  [9.16]  is  1.608  and  that  obtained 
by  using  the  exact  formula  [9.8]  is  1.600  , corresponding,  if 
A ^ = 1 , to  compressions  of  21.14  and  20.94  per  cent  respectively. 

A parallel  insensitivity  of  the  compressive  force  to  the 
precise  form  of  the  strain-energy  function  is  also  found.  The 
compressive  force  F , per  unit  cross-sectional  area  measured 
in  the  undeformed  state,  is  given  by 

F = -4A‘!5T0(a1n2+a2n1*  + bn6+cn8)  + 0(n10)  , [9.20] 


where 


and 


» ' J (Vif) 


. HL  (A  Hi) 

405  ^0  595;  ’ 


T 

0 


(W  +A^W  ) 
1 3 2 


A = 1 


[9.21] 


[9.22] 
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